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ABSTRACT 



Results of simultaneous INTEGRAL and optical ob- 
servations of galactic microquasar SS433 in May 2003 
are presented. The analysis of the X-ray and optical 
eclipse duration and hard X-ray spectra obtained by 
INTEGRAL together with optical spectroscopy ob- 
tained on the 6-m telescope allows us to construct 
a model of SS433 as a massive X-ray binary. X- 
ray eclipse in hard X-rays has a depth of ~ 80% 
and extended wings. The optical spectroscopy al- 
lows us to identify the optical companion as a A5-A7 
supergiant and to measure its radial velocity semi- 
amplitude K v = 132 km/s. A strong heating effect 
in the optical star atmosphere is discovered spectro- 
scopically. The observed broadband X-ray spectrum 
2-100 keV can be described by emission from opti- 
cally thin thermal plasma with kT ~ 15 — 20keV 

Key words: INTEGRAL; SS433; X-rays; optical 
spectroscopy. 



1. INTRODUCTION 



SS433 is a supercritically accreting microquasar with 
precessing accretion disk and collimated mildly rela- 
tivistic (v « 0.26c ) jets (see e.g. Cherepashchuk et 
al. 1996 and references therein). Since the discovery 
in 1978, this unique massive X-ray binary system 
has been investigated in optical, radio and X-rays 
by many authors (e.g. Cherepashchuk 1981, Mar- 
gon 1984, Kawai et al. 1989, Kotani et al. 1996, 
Goranskii et al. 1998a,b, Marshall et al 2002, Fab- 
rika 2004). 

SS433 is highly variable system and shows different 
types of periodicities: 



1. Precessional variability (P prec = 162 d .5) which is 
observed by periodical Doppler-shifts of H, He I, Fe 
XXV and other X-ray emission lines in optical and 
X-ray spectra and is clearly visible in optical and 
X-ray light curves. 

2. Orbital periodicity (P orb = 13 d .082). The shape 
of the orbital optical light curve strongly changes 
with precession phase (Goranskii et al. 1998 a,b, 
Cherepashchuk and Yarikov 1991). The value of 
the orbital period remains stable over more than 
30 years, which can be considered as an argument 
against the common envelope model for SS433. 

3. Nutation periodicity {P nu t — 6 d .2877) which is 
observed as periodic deviations from purely preces- 
sional Doppler motion of emission lines and can be 
recovered from photometric data (Goranskii et al. 
1988 a,b, Cherepashchuk and Yarikov 1991). The pe- 
riod and phase of the nutational variability remain 
stable over at least 16 years (around 950 nutation 
periods). Nutation radial velocity variations are de- 
layed from the nutation photometric variability by 
0.6 days, which corresponds to the travel time from 
the accretion disk center to the formation region of 
optical moving emission lines downstream the rela- 
tivistic jets located at I ~ 10 14 — 10 15 cm away from 
the disk center. 

The AO-1 INTEGRAL observations of SS433 discov- 
ered hard (up to 100 keV) X-ray spectrum in this su- 
percritically accreting microquasar (Cherepashchuk 
et al. 2003), suggesting the presence of an extended 
hot (with temperature up to 10 8 K) region at the 
central parts of the accretion disk. These new data 
made it possible to compare the eclipse characteris- 
tics of SS433 at different energies: soft X-rays (2-10 
keV, the ASCA data), soft and medium X-ray (1-35 
keV, the Ginga data), hard X-ray (20-70 keV, the 
INTEGRAL data), and optical. Such a comparison 
allows us to study the innermost structure of the su- 
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Figure 1. IBIS/ISGRI spectrum of SS433 collected 
over INTEGRAL orbits 67-69 in May 2003. The best 
fit (solid line) is for exponential cut-off with kT = 
14 ±2 keV (reduced x 2 = 0.32 for 7 dof, CL w 95%;. 
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Figure 3. IBIS/ISGRI 20-40 keV (upper panel) and 
40-70 keV (bottom panel) count rates of SS433 (with- 
out background subtraction). The egress part of the 
X-ray eclipse was kindly provided by Diana Han- 
nikainen. 




Figure 2. RXTE PCA and IBIS/ISGRI spec- 
trum of SS433 obtained during simultaneous 
RXTE/INTEGRAL observations of SS433 in March 
2004 (JD 2453076-2453078). The best fit (solid 
line ) is for emission from optically thin plasma with 
kT = 19.5 ± IkeV. 



percritical accretion disk and to constrain the basic 
parameters of the binary system. This is especially 
important in view of recent controversial results on 
the optical spectroscopy of SS433 (Gies et al. 2002, 
Charles et al. 2004). 



2. Crimean Station of Sternberg Astronomical In- 
stitute (V. Lyuty, T. Irsmambetova). Made BVR- 
photometry at the 0.6-m telescope. 

3. Kazan State University. Performed optical 
photometry at the Russian- Turkish 1.5-m telescope 
RTT-150 at the TUBITAK National Observatory, 
Turkey (I. Bikmaev, N. Sakhibullin) . 

4. Infrared K-photometry was performed at 1.1-m 
AZT-24 telescope of Pulkovo Observatory in Italy, 
Observatore di Campo Imperatore (Yu. Gnedin, A. 
Arkharov, V. Larionov). 

5. Radio monitoring at cm wavelength has been per- 
formed at the RATAN-600 radio telescope of SAO 
RAS, indicating the source to be at its non- flaring 
state (see Trushkin 2003 for more detail). 

The results of a quick-look analysis of the INTE- 
GRAL observations of SS433 were published earlier 
(Cherepashchuk et al. 2003). Here we present the re- 
sults of more detailed study of the INTEGRAL data 
together with the simultaneous optical spectroscopic 
and photometric observations of SS433. More com- 
plete studies of all data on SS433 obtained during 
this campaign will be reported elsewhere. 



2. OBSERVATIONAL CAMPAIGN 



The coordinated multiwavelength observational cam- 
paign of SS433 was organized during the AO-1 IN- 
TEGRAL observations of the SS433 field in March- 
June 2003 with the participation of the following re- 
search teams. 

1. Special Astrophysical Observatory of Russian 
Academy of Science (SAO RAS) (S. Fabrika, E. Bar- 
sukova). Provided high signal-to- noise optical spec- 
troscopy of SS433 at the 6-m telescope. 



3. IBIS/ISGRI LIGHT CURVES AND SPECTRA 

The INTEGRAL IBIS/ISGRI data were analyzed us- 
ing publically available ISDC software (OSA-3 ver- 
sion). The main results are as follows. 

3.1. The IBIS/ISGRI spectrum of SS433 (25-100 
keV) obtained in May 2003 can be equally well fit- 
ted both by a power-law ~ E~ a with photon in- 
dex a Ri 2.7±0.13 or optically thin thermal plasma 
emission with kT ~ 12 - 16 keV (Fig. Tlie 
source was not significantly detected by the Jem- 
X telescope during these observations. However, we 
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Figure 4- Upper panel: Averaged over 20 ks 
IBIS/ISGRI 25-50 keV eclipse light curve superim- 
posed on the mean V-light curve of SS433 (filled tri- 
angles) simultaneously obtained at the SAI Crimean 
Laboratory 0.6-m telescope, normalized to the max- 
imum flux. Bottom panel: The same IBIS/ISGRI 
light curve and normalized to the maximum flux 
Ginga 5-28 keV light curve (filled sguares; from 
Kawai et al. 1989) and ASCA 6-10 KeV light curve 
(filled triangles; from Kotani et al. 1996). 



formed simultaneously with INTEGRAL in the end 
of March 2004 (PI: M.Revnivtsev) to obtain the 2- 
100 keV spectrum of SS433. The source was observed 
about one binary orbit (~ 14 days) later than the 
disk maximum opening phase. The resulting spec- 
trum is shown in Fig. and is fitted by emission 
from optically thin thermal plasma with kT ~ 20 
keV. 



The integrated hard X-ray luminosity is La, (18 — 
60keV) ~ 4x 10 35 erg/s, 2^(60- 120keV) ~2x 10 35 
erg/s (assuming the 5 kpc distance to SS433), which 
is about 10% of the soft X-ray jet luminosity. 



3.2. No rapid variability was found in the hard X-ray 
band on characteristic timescales < 1 hour. 



3.3. The IBIS/ISGRI count rates of SS433 are pre- 
sented in Fig. |21 The X-ray eclipse at hard energies 
is observed to be a little narrower than the optical 
one, slightly broader than in the 1-35 keV energy 
range and display extended wings (Fig. 0J. This is 
opposite to what is found in ordinary eclipsing X-ray 
binaries (like Cen X-3, Vela X-l etc.), in which the 
X-ray eclipse duration decreases with energy. This 
new fact may reflect a complicated structure of the 
innermost supecritical accretion disk in SS433. 



3.4. The eclipse depth is about 80% in hard X-rays 
compared to ~ 50% in 6-10 and 5-28 keV band (Fig. 
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Figure 5. Best fit of the IBIS/ISGRI X-ray eclipse 
light curve by the binary system model with the mass 
ratio q = m x /m v = 0.3 (the solid curve) and q=0.75 
(the dashed curve). Hard X-ray emission is modeled 
by a semi- ellipsoid with axes dj — 0.3a, bj = 0.75a 
posed at the center of a thick disk with radius r^ — 
0.47a (for the case q = 0.75,) and aj — 0.3a, bj — 
0.75a, r d = 0.37a (for q = 0.3). Here a w 85R G (for 
q = 0.3/ is the binary semi-major axis. 



3.5. The 25-50 keV X-ray flux increases from ~ 5 
to ~ 17 mCrab during precession period in Marh- 
May 2003. This modulation is 1.5-2 times larger 
than observed in 2-10 keV egergy band. Thus, both 
precessional and eclipsing hard X-ray variabilities in 
SS433 exceed by 1.5-2 times those in the standard X- 
ray band. This suggests a more compact hard X-ray 
emitting region in the central parts of the accretion 
disk. 



The observed X-ray eclipse was interpreted by the 
model of a close binary with thick precessing ac- 
cretion disk (Antokhina et al. 1992). The model 
includes the following basic elements: i) Dark opti- 
cal star filling its Roche lobe, ii) Dark thick pre- 
cessing accretion disk with a central cone-like re- 
gion (the cone angle ~ 60°). iii) Bright extended 
bulge located at the central cone of the accretion 
disk. The observed X-ray eclipse in the 25-50 keV 
band is best fitted by this model for the mass ratio 
q = m x /m v ~ 0.2 and the bright bulge size ~ 0.3a, 
where a ~ 85_R© (for q = 0.3) is the binary semi- 
major axis (Fig. [5j. However, for reasons discussed 
below we shall assume q — 0.3 as a fiducial value. 
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Figure 6. V-light curve of SS 433 obtained at 
the RTT150 telescope (TUBITAK National Obser- 
vatory, Turkey) simultaneously with INTEGRAL ob- 
servations. In the bottom: photometry of control 
stars (V N3 = 12.975 ± 0.006, Vna = 12.746 ± 0.007,). 



4. OPTICAL OBSERVATIONS 
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Figure 7. JHK-photometrical light curve of SS433 
obtained by AZT-24 1.1-m IR telescope in July- 
August 2003 (Campo Imperatore, Italy). 



4.1. Optical photometry 

Photometry of SS433 was performed simultane- 
ously with INTEGRAL observations by the Russian- 
Turkish RTT-150 telescope of Kazan University at 
the TUBITAK National Observatory Turkey. Ob- 
servations were made by using commercial, Ther- 
moelectrically cooled to -60 C ANDOR firm CCD 
(model DW436, www.andor-tech.com) provided for 
RTT150 by MPA. Low noise back-illuminated EEV 
chip has 2048 x 2048 pixels of 13.5 mkm size each. 
Full field of view is 8 x 8 arcmin with frame read- 
ing time of 40 sec at 2 x 2 binning. To increase the 
time resolution only parts of the field of view with 
reference stars N 1,2,3,4,5 ( Leibowits and Mendel- 
son, 1982) around SS433 were stored to PC. The 
obtained V-light curve is presented in Fig. El Strong 
(~ 0.15 mag) intranight variability of the source on 
timescales ~ 100 s - 100 min is clearly detected. 
The mean V light curve during the eclipse obtained 
at the Crimean Laboratory of SAI simultaneously 
with the INTEGRAL observations is also shown in 
Fig. 01 The optical eclipse minimum is observed 
at JD = 2452770.863, as predicted by the orbital 
ephemeris given by Goranskii et al. 1998a. 



4.2. IR-photometry 

Near IR observations of SS433 were obtained at 
the AZT-24 1.1m telescope in Campo Imperatore 
(Italy) with SWIRCAM during a period spanning 



July-August 2003. SWIRCAM is the infrared cam- 
era that incorporates a 256x256 HgCdTe NIGMOS 
3-class (PICNIC) detector at the focus of AZT-24. 
It yields a scale of l".04/pixel resulting in a field of 
view of 4x4 sq. arcmin. The observations were 
performed through standard JHK broadband filters. 
The NIR monitoring of SS433 started several orbits 
after the INTEGRAL observations. The JHK light 
curves obtained in July- August 2003 are presented 
in Fig. [7| The observations were done close to the 
cross-over precessional phase, where the orbital mod- 
ulation appears to be significantly reduced. 



4.3. Optical spectroscopy 

The optical spectroscopy of SS433 was carried out 
at the 6-m telescope of SAO RAS during 6 nights in 
April 28 and May 9-13 2003, simultaneously with the 
INTEGRAL observations, at the precessional phase 
corresponding to the maximum disk opening angle. 
The long-slit spectrograph and PM-CCD 1000x1000 
detector were used to obtain 4100-5300 A spectra 
with a resolution of 3 A (1.2 A/pix). The blue part 
of the optical spectrum including He II 4686 emis- 
sion line was chosen for the analysis. The standard 
technique was used for spectral reduction and cal- 
ibration. Fig. |S1 shows four spectra of supergiants 
with known temperatures (Lc Borgnc ct al. 2003) 
and SS433 spectra averaged over nights 28/04/03 
and 11/05/03 (the disk eclipse phase) (second from 
the top spectrum in the figure) and over 09/05/03 
and 12/05/03 (corresponding to evenly spaced from 
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Figure 8. Optical spectra of SS433 taken on the 6-m telescope and spectra of supergiants of known tempera- 
tures (Le Borgne et al. 2003) for a comparison. Two upper spectra are averaged over nights 28/04/03 and 
11/05/03 (the disk eclipse phase) (second from the top spectrum in the figure) and over 09/05/03 and 12/05/03 
(corresponding to evenly spaced from the eclipse center orbital phases 0.88 and 0.12) (the top spectrum). 




Figure 9. Evolution of absorption lines in the optical spectra of SS433 taken at different orbital phases. For 
comparison, spectra of 8400 K and 20000 K supergiants are shown in the bottom. Orbital evolution of blends 
and radial velocities are traced by different absorption lines (dashes). See the text for more detail. 
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Figure 10. Radial velocity curve of the optical com- 
panion of SS433 obtained from 4 individual absorp- 
tion lines over 6 nights. 

the eclipse center orbital phases 0.88 and 0.12) (the 
top spectrum). Spectral resolution in SS433 spec- 
tra and in those of standard supergiants is the same. 
The signal-to-noise ratio in our spectra is > 60 at 
A4250 A per resolution element. Fe II 4233 emission 
is seen, emission lines in SS433 spectrum are broad 
{FWHM ~ 10 A). CII 4267 and NIII 4196,4200 
+ Hell 4200 emissions are also marginally present. 
In the Figure, the absorption lines seen out of the 
eclipse are marked with short vertical bars; those 
which are seen in the eclipse or present in both spec- 
tra are marked with long bars. 

The intensity of absorption lines during the eclipse 
allows us to estimate the effective temperature of 
the optical star to be T < 9000 K. The relative 
intensities of the strongest absorption lines indicate 
T = 8000 ±500 K, implying the optical spectral class 
of the companion A5-A7I. 

The heating effect of the companion atmosphere 
is discovered by absorption lines. During the 
disk eclipse egress, low-excitation absorption lines 
strongly weaken. The stellar hemisphere illuminated 
by the bright accretion disk probably has a temper- 
ature of ~ 20, 000 K, as the presence of CII 4267 
absorption suggests. This is the only strong line 
in this region in hot stars. The CII 4267 absorp- 
tion is deep out of the eclipse (Fig. |HJ and the 
line is only marginally detected inside the eclipse 
phases. Evolution of the blend 4273 (Fel 4271.8 + 
Fell 4273.3 + CrII 4275.6) is seen: in the eclipse, low- 
excitation line Fel is strongest, while out of eclipse 
Fell+Crll lines enhance. Also, other Fel lines ap- 
pear only in the eclipse. These effects are illustrated 
better by Fig. [5J in which spectra of standard stars 
with different effective temperatures are shown to- 
gether with spectra of SS433 obtained in 9/05, 10/05, 
28/04+11/05, 12/05 and 13/05 (orbital eclipses fell 
in nights 10/05 and 28/04+11/05). 

In Fig. 1101 we present individual absorption line ra- 
dial velocities for four best-suited lines averaged over 
one night, as a function of the orbital phase for 6 
nights. It is important to note that the strongest 



-200 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 

0. S 0. 8 1 1.2 1.4 

Orbital Phase 

Figure 11. Mean radial velocity curve of the optical 
companion of SS433 measured from 22 individual ab- 
sorption lines. The accretion disk radial velocities as 
measured from by Hell 4686 emission (Fabrika and 
Bychkova, 1990) is shown by the dotted line. 



lines (Fell, Hel, hydrogen) are present in SS433 spec- 
trum as emissions. Some other strong lines (among 
Fell and Hel lines) appear in SS433 as P Cyg-like ab- 
sorption components or distorted absorptions, they 
are disfavored for the optical star radial velocity anal- 
ysis, as they most probably form in the powerful out- 
flowing disk wind (Fabrika et al. 1997). 

The radial velocity curve measured by the most reli- 
able collection of 22 absorption lines in the spectral 
range 4200-5300 A is shown in Fig. EI Note that the 
absorption lines in SS433 spectrum are seen at our 
spectral resolution as blends containing 1-3 strong 
lines. Relative intensities of the lines are probably 
change from date to date through the eclipse. In 
this reason we measured only those lines and in only 
those dates where the lines were detected the most 
convincingly. The radial velocity (Fig. Illfl has been 
obtained by coadding of radial velocity curves of in- 
dividual lines and a radial velocity of all lines for 
11/05/03 (the orbital phase 0.048) was adopted to be 
a zero. The derived radial velocity semi-amplitude 
of the optical star is K v — 132 ± 9 km/s, the 
gamma- velocity of the binary system is v 1 = 14 ± 2 
km/s. The absorption line radial velocity transition 
through the 7-velocity occurs at the middle of the 
optical eclipse (<j>b = 0.07), confirming that the lines 
actually belong to the optical star. 

These results confirm the earlier determination of K v 
by Gies et al. (2002) also made at the maximum disk 
opening angles (note that spectroscopic observations 
by Charles et al. 2004 were performed at the cross- 
over phase of SS433 when the accretion disk is seen 
edge-on; such a phase is disfavored for the optical 
star radial velocity analysis as strong gas outflows are 
present in the disk orbital plane; selective absorption 
in this moving gas affects the true radial velocity of 
the optical star). 

The heating effect of the optical star also distorts the 
radial velocity semi-amplitude. The analysis (Wade 



that true value of the semi-amplitude of the radial 
velocity curve as derived from these absorption lines 
can be reduced to ~ 70 — 80 km/s. 



5. BASIC PARAMETERS OF THE BINARY 
SYSTEM 



Comparison of radial velocities of the accretion disk 
(K x = 175 km/s, Fabrika and Bychkova 1990) and 
optical star (K v = 132 km/s) yields the mass ratio 
in the SS433 system q = m x /m v — K v /K x — 0.75. 
This is an upper limit, as the true value of K v 
should be corrected for the heating effect. Taken at 
face value, these q, K v lead to the optical star mass 
function f v w 3.12M and the binary component 
masses m x — 18M Q , m v — 24A/ Q . However, such a 
big mass ratio is in a strong disagreement with the 
observed duration of X-ray eclipse, which suggests 
much smaller mass ratio q ~ 0.2 — 0.3. We stress that 
the binary inclination angle in SS433 (i = 78°. 8) is 
fixed from the analysis of moving emission lines. 

There are two possibilities: (1) either the model we 
used to fit X-ray eclipses should be modified, or (2) 
the value of K x and K v are influenced by additional 
physical effects. Though there are some reasons to 
modify the model (e.g., asymmetric shape of the X- 
ray eclipse, which may suggest an asymmetric wind 
outflow from the illuminated part of the optical star) , 
we shall consider here only the second possibility. We 
repeat again that the actual value of K v should be 
decreased to account for the observed heating effect. 

Let us assume the mass ratio in the system to be 
q = 0.3, as with this value we can satisfactorily de- 
scribe the X-ray eclipse width (see Fig. 03 . Taking 
K v = 132 km/s yields /„ « 3.1M Q and m x = 62M Q , 
m v = 2O6M , K x = 440 km/s. Clearly, this is an 
unacceptable model. Now let us decrease K v down 
to 85 km/s, the lower limit which follows from more 
accurate treatment of the heating effect in the radial 
velocity curve analysis (Wade and Home 1988, An- 
tokhina et al. 2004). This would yield a better fit 
with m x — 17Mq, m v = 55M Q , and K x = 283 km/s, 
still too high to be acceptable. Finally, let us take 
into account that less than half of the stellar surface 
is sufficiently cool to give the absorption lines under 
study (e.g. due to sideway heating from X-ray jets 
and scattered UV radiation in the strong accretion 
disk wind). This additionally decrease the value of 
the actual radial velocity semi-amplitude. So tak- 
ing, for example, K v — 70 km/s and q = 0.3 yields 
the optical star mass function /„ = 0.46M Q and bi- 
nary masses m x = 9M , m v — 30M© and optical 
star radius R v ~ 4Oi? . This radius is compatible 
with typical bolometric luminosity of a 30M Q A5- 
A7 supergiant with T e ff ~ 8500 K. In this solution, 
K x ~ 233 km/s, larger than the measured value 175 
km/s, but the true value of K x may be affected by 
the strong accretion disk wind. 



The INTEGRAL and optical observations of SS433 
provide support that SS433 is a massive black hole 
X-ray binary system. The presence of broad wings 
in the observed hard X-ray eclipse profile evidences 
for an extended central hot region in the SS433 ac- 
cretion disk. This region is observed to be substan- 
tially eclipsed during precessional motion of the ac- 
cretion disk. The detected broad-band (2-100 keV) 
X-ray spectrum is best-fitteed by thermal plasma 
with kT ~ 15 — 20 keV. The origin of such an ex- 
tended hot region in the cenral parts of supercriti- 
cally accretion disk in SS433 needs more studies. 

The measured semi-amplitude of the radial velocity 
of the optical star is K v = 132 ± 9 km/s, which con- 
firms the results by Gies et al. (2002). However, the 
actual semi-amplitude should be decreased down to 
~ 70 km/s to be in agreement with the observed du- 
ration of X-ray and optical eclipses. This correction 
seems indeed inevitable because of strong heating ef- 
fect present in the optical star of SS433. 

To measure more precisely the binary parameters 
of SS433, further hard X-ray observations of SS433 
are required at different precession phases and high- 
resolution optical spectroscopy need to be carried out 
at the maximum disk open angle to accurately mea- 
sure K v and K x . Analysis of X-ray observations at 
minimum X-ray flux and the cross-over precession 
phase will allow more detailed reconstruction of X- 
ray emitting zone in the center of the supercritically 
accretion disk in SS433. 
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